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t""^ ■ We present solutions describing spatially closed, open, or flat cosmologies in the 

massive gravity theory within the recently proposed tetrad formulation. We find that 
the effect of the graviton mass is equivalent to introducing to the Einstein equations 
g_) ' a matter source that can consist of several different matter types - a cosmological 



term, quintessence, gas of cosmic strings, and non-relativistic cold matter. 
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^ . The currently observed acceleration of our universe [l| is the main motivation of attempts 



to try to modify the theory of gravity, for example by giving a tiny mass to the graviton. 
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t^ ■ This can effectively give rise to a small cosmological term, at least within the simplest 

bimetric models of massive gravity j^.^The theory of massive gravity is not unique and 



there exist a number of its models (see (3| for a review). Recently a massive gravity model 



^ . was proposed which could be in some sense special, since it is the only one that admits a 

pi ! quadratic action in the unbroken phase [^. This model reproduces in fact that in 5| but 

uses a parametrization which avoids dealing with the square root of a tensor, so that it is 
much simpler to work with. In what follows we use the formulation of jj| to construct new 
cosmological solutions, in particular those with spatially open and closed metrics. It turns 
out that the graviton mass gives rise not only to a cosmological term, but can also manifest 
itself as other matter types, with different equations of state. 

The new parametrization of the massive gravity uses the tensor jj| 

Sab = e'Xd^cpB - Vab , (1) 

where 0b are four scalars, tjab is the Minkowski metric, the tetrad e^ determines the metric 



gij.iy — ^AB^i^^i^ g^^j jg constrained to fulfill the conditions 

eAd^(j)B = e'^d^cpA , (2) 

which insure that 5*^^ is symmetric. The action can be chosen to be 

1 /■ / 1 



S = ^J i-^R + m^CJV^d^x + Srn, (3) 

here Sm describes the ordinary matter (for example perfect fluid), m is the graviton mass, 

r ^{q2 qAqB\ I ^3 ABCQqMqNqP I ^4 ABCD qM qN qP qQ f a\ 

^ - 2^^^ ~ ^^B^a) + -^(^MNPQ^ ^A^B^C + J^^MNPQf^ bAO^b^^jj, [4) 

where 5" = Sa and C3, C4 are two parameters (for c^ = C4 = the action becomes quadratic). 
This theory is equivalent to the one proposed in (5|, where it was parameterized by Kj^ = 
b^ + \/ d^<\)Adv<\>^- Since Sab = — e^e^i^^j,, the two parameterizations are equivalent, but 
the Sab parameterization has the obvious advantage, since it does not require to take the 
square root which is only defined through a series expansion. 
The field equations of the above action are 



G,. = m%, + SttGT^-, (5) 



with 
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T^, = {S + l)ep,<PA - d,<PAdA'' + ^ eABCQe'''''''^e^^dAM S^Sp 



+ J eABCDe'''^'''^e^d,^M S§S^Sg - g,X , (6) 

which is symmetric in view of ([2]) and which has to satisfy the conservation condition 

V^T^. = 0, (7) 

also V^T™ = 0. Let t, r, t?, </) be spherical coordinates and n'^ = (sin ^ cos (f, sin t? sin (f, cos •&) 
a unit vector. We choose the four scalars to be 

0° = T(t,r), r = U{t,r)n\ (8) 

and the tetrad 

ejrfx^ = g(t, r)dt + h{t, r)dr, e^dx^" = ^f^d{p{t, r)n''). (9) 

p[t, r) 



Choosing h = E?pp' / {qp^) the metric becomes diagonal 



S)f.^--^" 



ds'=[l-^-^\[ q^df - -UL. dr^ ) - R\d^' + sin^ Mcp'). (10) 



Inserting ([8]), ([9]) into ([2]) , ([5]) , ([7]) , the angular variables decouple, giving a system of non- linear 

PDE's for T(t, r), U{t, r), q{t, r),p{t, r), R{t, r). 
Let us consider first the static case. Setting 

/ f N(r) \ 
qit,r) = q{r), R(t,r) = r, p(t, r)=exp( / dr \ , (11) 

the metric becomes 

ds^ = q^de - NHr^ - r^ {dd^ + sin^ M^^). (12) 

If T™ = then the field equations are fulfilled by 

f-^.-^-'-^-'^^r^ (13) 

and T(t,r) = t, U{t,r) = Cr, provided that 

2-C 3-3C + C2 

C3 = ^r^I' ^4 = - (c-i)2 > (14) 

where C, M are integration constants. This solution was in fact found in [6| , it describes the 
Schwarzschild-dS or (Schwarzschild-AdS) metrics. 

Let us now consider time-dependent solutions. Choosing 

g(t,r) = a(t), Rit,r) = ait)r, p{t,r) = Z' , (15) 

1 + vl — i^f 



with iiT = 0, ±1, the metric becomes 



dr 
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ds' = a'(t) ( de ^-— - r\d§^ + sin^ ^rfc^^^ . (16) 

\ 1 — Kr^ J 

The matter source is chosen to be SttG'T™';^ = diag(p(t), —P{t), —P{t), —P{t)), whose con- 
servation condition is 

p + 3-(p + P) = 0. (17) 

Choosing U{t,r) = Cra{t) and 

K = 0: T{t,r) = -^Cr^d 



K = ±1: T{t, r) = KCa\/\ - Kr'^ 

(18) 



we find that the conservation conditions (jTj) and the symmetry conditions ([2]) are fulfilled, 
provided that C3, C4 are again constrained according f lT^ . while Einstein equations ([5]) reduce 
to 

3^^±^ = m^(C-l) + p. (19) 

This describes the late time acceleration. Indeed, if p = 7P then p ~ a"^"^/''' so that for 
large a the second term on the right in (fT9|) becomes negligible. The dynamic is then driven 
by the first term corresponding to the cosmological constant im?{C — 1) which can be positive 
or negative, depending on value of the integration constant C . It is worth noting that the 
K = cosmology was previously described in the literature in the decoupling limit (5|. 
One can also get solutions for generic values of C3, C4. With f/(t, r) = Cr and 

K = {]: p{t,r)=r, T{t,r) = 



K = -1: p{t, r) = re\ T{t, r) = CVl + r 
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K = 1: pit, r) = re'\ T{t, r) = -iCVl^^ (20) 



the metric again becomes flT6|) . while the field equations reduce to 

a^ + Ka^ ,f ^ , 3C(3 + 3C3-C4) , 3C\c, - 2cs - I) , C^c^ ~ c,) , , 
3 7 = m C4 - 4c3 - 6 H 1 \ +p. 



a^ \ a a? a^ 

'(21) 

Although the matrix (9^0^ is degenerate for these solutions, this should not affect the dy- 
namics of their perturbations, so that the gravitons should still be massive. We see that the 
value of the cosmological term 772^(04 — 4c3 — 6) is now determined by the model parameters 
and not by an integration constant. Secondly, the graviton mass also gives rise to the three 
additional terms in the right hand side of the equation which decay as 1/a, 1/a^ and 1/a^. 
These can be viewed as three different matter types contributing to the energy density. The 
1/a^ behavior of the energy is the same as for the non-relativistic cold matter (P = 0), the 
1/a^ is the same as for a gas of cosmic strings, while the 1/a is typical for quintessence. 
We also notice that the K = 1 solution in (120|) . obtained via analytic continuation of the 
K = —1 case, has complex e^ and (pA, although the metric, S^b and the equation (l^Tj) are 
real. It seems that there should be a different representation of this solution, with real p{t, r) 
and T(t,r), but we could not find it. 

Summarizing, we have found spatially open, closed, and fiat cosmologies with massive 
gravitons. The graviton mass gives rise to a cosmological term thus leading to the late time 



acceleration, but it can also give additional contributions to the energy that can be viewed 
as matter components with different equations of state. 
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